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Abstract 

We consider the collider signature of right-handed neutrinos propagating in 5 (large) extra 
dimensions, and interacting with Standard Model fields only through a Yukawa coupling to the 
left-handed neutrino and the Higgs boson. These theories are attractive as they can explain the 
smallness of the neutrino mass, as has already been shown. We show that if 5 is bigger than two, it 
can result in an enhancement in the production rate of the Higgs boson, decaying either invisibly or 
to a b anti-6 quark pair, associated with an isolated high px charged lepton and missing transverse 
energy at future hadron colliders, such as the LHC. The enhancement is due to the large number 
of Kaluza-Klein neutrinos produced in the final state. The observation of the signal event would 
provide an opportunity to distinguish between the normal and inverted neutrino mass hierarchies, 
and to determine the absolute scale of neutrino masses by measuring the asymmetry of the observed 
event numbers in the electron and muon channels. 
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I. INTRODUCTION 



The Standard Model (SM) of high energy physics suffers from the gauge hierarchy 
problem, which is the fine tuning required to maintain a low electroweak scale (Mew ~ 
10 3 GeV) in the presence of another seemingly fundamental scale, the Planck scale (the 
scale of gravity, M v \ ~ 10 19 GeV). Supersymmetry, technicolor and more recently extra 
(space) dimensions have been proposed to address the hierarchy problem. 

Recent neutrino oscillation experiments have suggested a nonzero neutrino mass, with the 



best fit values of the mass differences and mixing angles given by 

A™Lar = 7 x 1(T 5 eV 2 , tan 2 9 solar = 0.4 



DDDD 



3 „\t2 +™2 , 



Am 2 tm = 2.5 x l(T d eV 2 , tan^ otm = 1 . (1) 

The oscillation between the three active flavors, v^v^v T ^ in the SM, accommodates this 
satisfactorily, with Am 2 olar = (m\ — m 2 ) and Am 2 tm = \m\ — m\\, where the are the 
physical neutrino masses. If the mj are also assumed to be of the same order of magnitude 
as the mass differences, it is quite challenging to explain why it is that the neutrinos are so 
light compared to the other leptons. 



It has been shown jy, yj that if there are other Large Extra Dimensions (LED) in addition 
to our usual four space-time dimensions, we could potentially solve the gauge hierarchy 
problem. It was then pointed out 0, that the smallness of the neutrino mass is 

naturally explained^ if right-handed neutrinos that propagate in some 5 number of these 
extra dimensions are introduced. We will refer to such neutrinos, which are SM gauge 
singlets, as "bulk neutrinos", as is the usual practice. Various aspects of theories with bulk 



neutrinos have been analyzed in Ref. 



Hi- 



We will take the view, as in Ref. that the standard three-active-flavor oscillation 
explains the data in Eq. (|T|). and that the mixing to sterile bulk neutrinos are small enough 

* In this work, we will not address the Liquid Scintillator Neutrino Detector (LSND) result Q. 
t We note here that the conventional see-saw mechanism to explain the smallness of the neutrino mass is 
equally appealing, but we will not consider it in this work. 
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to evade experimental constraints. In our earlier work Q|, we considered the constraints 
on theories with bulk neutrinos coming from neutrino oscillation experiments and also from 
requiring that perturbative unitarity be maintained in the theory in Higgs-Higgs scattering. 
We showed that strong constraints result when 5 > 2, though in that case, a precise 
calculation of the bound was not possible owing to the sensitivity on the cutoff scale, implying 
a dependence on the completion of the extra dimensional (effective) theory that we work 



12j , an alternative approach |8|, |10( wherein 



with. There, we noted that as pointed out in Ref. 

Eq. (P) is explained by the oscillation of the active species predominantly into sterile bulk 



neutrinos, appears to be disfavored by the Sudbury Neutrino Observatory (SNO) 
current data. 



2] 



neutral 



If bulk right-handed neutrinos are responsible for the smallness of the neutrino mass, we 
ask in this paper, what the consequences at a collider might be. We will show that in certain 
cases, the Higgs production and decay could be measurably altered from the SM expectation. 
We will present the discovery potential of hadron colliders, mainly focusing on the CERN 
Large Hadron Collider (LHC), through the mode: qq' — > W — > Higgs + lepton+ Br • The 
invisible decay mode of the Higgs boson in theories with a bulk right-handed neutrino at 



hadron colliders have also been explored in Refs. 



14 



The rest of the paper is organized as follows. We will introduce the extra dimensional 
theory with bulk neutrinos in Sec. |Hj write down the equivalent four dimensional Kaluza- 
Klein theory, with particular focus on the interaction of the right-handed bulk neutrino with 
the Higgs field and the left-handed neutrino. In Sec. IIII| we will analyze the collider signature 
of such a theory at the Fermilab Tevatron and the CERN LHC by performing a Monte-Carlo 
simulation, and improve the significance by appropriate cuts. Our conclusions are given in 
Sec.HVl 



II. RIGHT-HANDED NEUTRINOS IN EXTRA DIMENSIONS 

In this section, we summarize the framework with bulk right-handed neutrinos, details of 
which are given in Ref. jl3| . 
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To address the gauge hierarchy problem, Arkani-Hamed, Dimopoulos and Dvali (ADD) [6] 
postulate that the Standard Model (SM) fields are confined to a 4 dimensional (4-D) sub- 
space (brane) in an extra dimensional world of 4 + n total dimensions. ADD take the view 
that the only fundamental scale in nature is M*, which is of the order of Mew-i an d the 
apparent 4-D gravity scale (M p i) is then given by 

Ml = Ml +n V n , (2) 

where, V n is the volume of the (compact) extra dimensional space. In the simple case of 
each of the compact extra dimensions being of equal radius R', we have V n ~ R' n . Thus 
ADD argue that M p i appears to be a large scale from a 4 dimensional perspective simply 
because the volume V n is large. In other words, the explanation of why M p i is large is recast 
to stabilizing R' at a large value, so that V n is large. 

It should be pointed out that for a given M*, if the n compact dimensions have equal 
radii R', Eq. (J2J implies a particular value of R'. However, if it happens that there are two 
sets of compact extra dimensions of unequal size, 5 of them (5 < n) with radius R, and the 
other (n — 5) with radius R', then we have in this case, V n ~ R' {n - 5) R 5 . We can in this case 
think of R as an independent variable with R' being determined by Eq. (J2J). 

We consider the ADD framework, to which is added three (one for each generation) bulk 
fermions, ^ a {x li ,y), that propagate in 4+5 dimensions (5 of them compact with radius R), 
where the indices a, (3 — (1, 2, 3) denote the three generations, and y stands for {y 1 , ...,y 5 }. 



A. The Lagrangian 



We can split the Lagrangian into a bulk piece and a brane piece, 



S = I d 4 x d 5 y £ Bu ik + 5(y)C Br 



(3) 



^Buik contains the Einstein-Hilbert bulk gravity term (which we will not show explicitly, but 
can be found, for example, in Ref. the kinetic energy term for the bulk neutrino field 

\l/(x M , y) and in general, a bulk Majorana mass term for which for simplicity we will omit 



(see Ref. |8j for implications of a nonzero bulk Majorana mass). ^Brane contains the SM 
Lagrangian plus an interaction term between SM fields and ipR, 

£ Bu ik D * a iT M D M y a , 

^Branc D ^SM ~ I ~ F=^= hl^tf + h ' C - I > ( 4 ) 



where, is an 0(1) Yukawa coupling constant. It should be kept in mind that ipR is a 
function of {x^,y) whereas the SM fields are functions of only. The index M runs over 
{x»,y}. 

We can perform a Kaluza-Klein (KK) expansion of the 4+5 dimensional theory and obtain 
an equivalent 4-D theory by writing, 

r R {x^y) = Y. VS^V) / e (k), (5) 

n 

where, n = (ni, ...,ns) is a vector in "number space", ip^ are the KK modes and fn(y) is a 
complete set over y. A similar expansion is made for t/>£. To reduce clutter, we will simply 
write n and y for n and y, respectively. We will use the notation n = (0, 1, ...) and h = (1, ...) 
(n excludes 0). f n is an orthonormal set, 

1R d 5 y f:(y)f m (y) = S nm i (6) 



o 



and a convenient choice is 



e r e r 



fn(y) = y<? = -7f7=, (7) 

with Ss-i the surface "area" of a unit sphere in 5 dimensions, and V$ = ^j^-R 6 is the volume 
of the extra dimensional space. 

We define the fields* Q, 

v r — Wr j 



* The other linear combinations {^gf^ — tp^ and ^V™^ — V'"^ are decoupled from the SM fields, 
and we will not consider them further. Also, ^ 
that it is excluded from the particle spectrum. 



and we will not consider them further. Also, with orbifold compactification, we can project out ifi^ so 



5 



a(n) 
R 



a(— n) N 
R I ' 



We substitute the KK expansion for the bulk fields ip^ and into Eq. 
equivalent 4-D theory, 



to get the 



a=l h 



I n I 



(u$%?» + h.c.) 
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-E 

a,/3=l 



m 



a/3 



h.c. 



where, |n| = ynf + ... + nj. We note here that in there is a tower of KK states 
(vj^\ v^) with Dirac masses approximately equal to \n\/R. With SU(2) broken by the 
Higgs mechanism, by the Higgs field acquiring a vacuum expectation value (VEV), (h) = v, 
we have the neutrino mass matrix given by, 



m 



Kpv 



S S - 



Kpv 



(9) 



and m^P can be much smaller than v ~ M^w-, if V$ is large for being 0(1). Henceforth, 
we will assume that unless noted otherwise, repeated generation indices a,/3,i,j, and KK 
indices n, h, m, are summed over. 



B. KK states 

For 5 > 1, the state with mass ^ can be degenerate, and we denote the degeneracy at 
the n th level by d n . (Strictly speaking, we should denote this as d\ n \, but we will just write 
this as dn.) For example, for 6 = 3, the state with mass 1/R has d\ = 3, corresponding to 
(■hi, 77,2, n 3 ) — > (1, 0, 0), (0, 1, 0) and (0, 0, 1), all of which have the same mass. For large \n\, 
the leading power dependence of d n in 5 extra dimensions is given by d n = Cftln^ -1 , where 
the c n are 0(1) numbers. We define do = 1. 

For large \n\, we can think of the hi as a continuum and the leading behavior is given by 
the surface of the (5 — l)-sphere of radius \h\ in number space, 

d n ~ Sg— i\h\ 6 ~ l (in 5 dimensions). (10) 
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4 =1 = 1, d 5 = 2 ~ 2vm, 4= 3 ~ 4vm 2 . 

For example, for 5 = 3, ~ 47r|n| 2 , which is the surface of the 2-sphere with radius \n\. We 
will often use the continuum approximation for estimating various quantities. 

In a collider, the heaviest KK state that could be produced in the final state is limited 
by the center-of-mass energy ^fs. We define N s to be the radius of the biggest sphere in 
{rii} space such that N s /R = ^/s. The sum over the KK states of certain quantities can be 
divergent and can depend on N s . We will elaborate more on this later, and we will see that 
the production rate of the KK states in association with a Higgs boson can depend on N s , 
especially strongly for 5 > 3. 



< = i m 4 , 


"S = 


(r a T4 , 




a(n) 


_ ai Mn) 
— ' V L 


a _ jod H 




(C)*4 , 


p a — 
e R — 



C. Higgs boson interaction 

We can make the Yukawa coupling in Eq. (JHJ) diagonal in generation space (a, (3) with the 
rotations 



a(h) 



where the unitary matrices / and r are chosen to diagonalize m^, so that, (r m )* m" 13 (0) = 
rriiS 1 ^ . Similarly, l e and r e are chosen to diagonalize the electron- type mass matrix. In the 
usual way, the charged current interactions now become proportional to the MNS matrix 

V MNS = l\l. (11) 

We can work in the basis in which the charged lepton mass matrix is diagonal. In this 
case we have from Eq. (JTTJ), Vmns — I- Based on the Vmns obtained from a global fit to 
the oscillation data, we take the / to be such that the solar neutrino mixing angle between 
v e & h>n is given by tan 2 6 efl = 0.4, and the atmospheric oscillation mixing between & v T 
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is maximal. § This implies 



/ 0.847 0.531 \ 

-0.376 0.599 0.707 
V 0.375 -0.599 0.707/ 



(12) 



We take the value of such that the lowest mass eigenvalue m, obtained by 
diagonalizing the neutrino mass matrix is consistent with Eq. (JTJ. Defining £j = rriiR, 
we showed in Ref. |13| that in order to satisfy the limits on the probability of an active 
neutrino oscillating into the sterile KK states, we are led to restrict ourselves to the case 
when En %d h < 1 (for all 5). 

We will show in Sec. IIII Al that the production rate of the KK neutrinos in association 
with the Higgs boson is proportional to the matrix 



m. 



= m Z^ 



(13) 



with mo the diagonal 3x3 physical neutrino mass matrix . We will find it convenient 
to work in a basis in which we write the Higgs interaction term in Eq. (jHJ) as 



m 



i/3 



+ h.c. 



(14) 



where we have rotated the z/^ with the matrix r while retaining the vl in the flavor basis, 
in order to explicitly keep only the physical matrix /, while absorbing the unphysical matrix 
r into the definition of u' R . In particular, we will also show in Sec. IIII Al that the total 
production rate (sum over e, [i and r channels) of our signal process is proportional to 



m 



Presently the combination of various neutrino oscillation experiments constrain |Am^ tm | 
leaving its sign undetermined, while the solar MSW effect j^] fixes the sign and magnitude 
of Ar?7.g 0lar . Furthermore, oscillation experiments are only sensitive to the mass differences 
due to which the absolute scale of neutrino masses is not known. This leads to more than 
one possibility for the values of the neutrino mass eigen- values, and in Ref. [3], we noted 



For simplicity, we assume the small mixing angle #13 ~ 0, which leads to I ~ 0. 



8 



m. 



m„ 



m, 



NORMAL 



m, 



m, 



INVERTED 



DEGENERATE 



FIG. 1: Neutrino mass schemes. 

the three neutrino mass schemes shown in Fig. ^ In each scheme, we give below the matrix 
m v which specifies the Higgs interaction that we use to calculate the collider signature in 
the next section. 



(i) Normal Hierarchy: mi « 0, m 2 ~ 0.008 eV and 777,3 ~ 0.05 eV, and 

/ \ 

m„ = 0.0045 0.005 -0.005 . 
V 0.036 0.036 j 

(ii) Inverted Hierarchy: mi, m 2 ~ 0.05 eV and m 3 « 0, and 

/ 0.043 -0.019 0.019 \ 



(15) 



777, 



0.027 

V 

(iii) Degenerate: mi,m 2 , m 3 all at some mass scale less than the limits discussed in Ref. 



0.03 




-0.03 




(16) 



Here, to illustrate the character of the signature, we take the three masses to be around 
1 eV, an arbitrary choice, and 

/ 0.847 -0.376 0.375 \ 



777, 



(17) 



0.531 0.599 -0.599 
V 0.708 0.708 J 

As shown in Eq. ()14j) . the Higgs boson couples to the tower of KK neutrino states with 
the same Yukawa coupling that is proportional to the small neutrino mass. When 1/R is 
small and 5 is large, there can be a large number of KK states that may overcome the small 
Yukawa coupling and generate an observable collider signature in the production and the 
decay of the Higgs boson. We shall discuss this possibility in the next section. 




FIG. 2: Signal process 

III. COLLIDER SIGNATURE 



A. Production of Higgs boson 



The right-handed neutrino does not carry any electroweak quantum numbers and therefore 
can be produced in a collider only through the Yukawa interaction. Although the production 
cross section is suppressed by the small Yukawa coupling (since it is proportional to the 
small neutrino mass), the cross section could get a large enhancement when summed over 
the large number of allowed KK excitations of the bulk neutrino. In such a case, the Higgs 
phenomenology can be altered in an interesting way and we study this in the process: 

qq'^W* -+£+hi/$ n \ (18) 

as shown in Fig. |21 where £ + is a charged lepton, v'^ 1 the n-th KK mass eigen-state and 
the process is mediated by z/£ with a — I. The production of v'r in the final state would 
lead to missing energy in the event. 

The total cross section for this process at a hadron collider is 

a{Pl P 2 ^ thvf) = $> (n) (PiP 2 -> thv$>) d n (19) 



dXxdx 2 
n q,q' 



■{x x <-> X 2 ) 



where J2n is a sum over the allowed KK states (0 up to N s = ^/sK) with deg eneracy d^i at 
the n th level, P\,Pi represent the hadronic initial state, f q /p(x,fx) is the parton distribution 
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function (PDF). CTEQ6 M 18] is adopted in our calculation. We take the factorization scale 
to be the invariant mass of the constituent process in our numerical calculation, is 
the parton level cross section to produce the n th KK mass eigen-state v'^ n \ and is given by 



a 



(<<) 



2 ydn 3 £X> 

2 spin 
color 



(n) 



where the bar above the 



M in) 



(20) 



denotes averaging over the initial state spin and color, dn 3 
represents 3-body final state phase space, and the squared matrix element resulting from 
Eq. (|14|) is given by 

2 



\MW\ 



1 1 

2 2 



in 



it 



1 



(p 2 - m 2 w 



2\2i 



m 



2\2 



x \ 32{pe+ ■ p q ) (p h ■ p ( n) ) (p h -Pq') + 32m V) (p e + ■ p q ) (p h ■ p?) 



(21) 



where p q and pq are the 4-momenta of the incoming partons q and q, pi+, ph and p ( n ) are 
the 4-momenta of the outgoing particles, p = p q + p^ is the momentum of the virtual W 
boson and q = Ph is the momentum of the intermediate z/£. The factor Q |J is due 

to the spin averaging, | is the color factor, and v(= 246 GeV) is the usual Higgs field VEV. 

In the continuum approximation (c.f. Eq. ()10|)) after summing over the allowed KK states, 
the total cross section (including e + ,/i + and r + channels) scales as 

Al\2 



\m. 



5' 



(22) 



Using Eqs. (fTU|) - (j2*T|) we compute the production cross section of the Higgs boson for different 
choices of m h , 1/R and 5 at the Fermilab Tevatron (pp collision at ^fs = 1.96 TeV) and the 
CERN LHC (pp collision at ^/s = 14 TeV). The sum of the Higgs boson production cross 



section (including 



is shown in Fig. El and Table HI for a representative 



choice of J2i.e \fn l y\ 2 = (0.05 eV) 2 . (The actual numbers for a particular mass matrix can be 



obtained easily using the scaling relation shown in Eq. (J22J).) Although we do not strictly 
impose the constraints on 1/R that we derived in our early work kJ, we choose the range 
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FIG. 3: Total cross section (with no cuts applied) as a function of 1/R for = 115GeV and 
m h = 150GeV at the Tevatron and the LHC for £^ |m*f | 2 = (0.05 eV) 2 . 



of keeping this in mind. For this range of 1/R, we see from Table H] that the signal rate 
would be too small to be seen at the Tevatron, and we therefore only consider the LHC in 
the rest of this work. Due to the summation over the KK states in Eq. (|19p. the production 
cross section increases with decreasing 1/R, owing to the increase in the number of KK 
states (c.f. Sec. Ill Bj) that can be produced as final state particles. Also, as the absolute 
neutrino masses m 8 increase, the ffi^ increase (for example see Eq. (|T7|) for the degenerate 
mass scheme) and the total cross section increases with the dependence shown in Eq. ([22)1 . 
However, as the absolute neutrino masses increase, the constraints on 1/R from neutrino 
oscillation data and from unitarity becomes stronger, and imposing the stronger constraint 
keeps the total cross section from increasing. 
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TABLE I: Cross section (in fb, with no cuts applied) of the process qq' — > l + hvn, with £ + = 
e + ,/i + ,r + , at the Tevatron and the LHC for different choices of nih, l/R and 5, for X^l^tl 2 = 
(0.05 eV) 2 . 





l/R 
(eV) 


m h = 115 GeV 


m h = 150 GeV 


Tevatron 


LHC 


Tevatron 


LHC 


5 = 3 


10 2 
10 3 
10 4 


49.2 
48.2 x 10~ 3 
49.2 x 10~ 6 


4.03 x 10 3 

4.03 
4.03 x 10" 3 


21.9 
21.9 x 10~ 3 
21.9 x 10~ 6 


2.92 x 10 3 

2.92 
2.92 x 10~ 3 


5 = 2 


1 


0.23 x 10" 3 


0.009 


0.093 x 10~ 3 


0.005 



In what we call the degenerate mass scheme, we arbitrarily pick the neutrino mass scale 
to be around 1 eV. Using the scaling relation in Eq. (J22J), one can easily get the total cross 
section for the degenerate mass scheme from Fig. |3] and Table HI In the following, when we 
present results for the degenerate case, we keep in mind the stronger constraints on l/R (see 

n 

Ref. [13]), and choose l/R such that we obtain the same order of magnitude cross section as 
for the normal and inverted mass schemes. 

We will study the decay modes of the Higgs boson in order to identify the collider signature 
of the signal events. In addition to the decay modes present in the SM, the Higgs boson can 
also decay into (is L + v^), an invisible decay mode. We will study this next. 



B. Decay of Higgs boson 

The invisible decay width (summed over all neutrino flavors) of the Higgs boson (h 
vl^r) is given by 



Tinvis = ^ndn = TZ 



- if 

m,. 



rn 



m h 1 



in) 



(23) 



where the sum over the KK states is up to N m = m^R since the most massive KK state that 
the Higgs boson could decay into is limited by m^. Using the continuum approximation to 
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TABLE II: Decay branching ratio of the Higgs boson for different choices of m^, 1/R and 6, for 
E M |mjf| 2 = (0.05 eV) 2 . 







1 / R 


m h - 


- 1 1 tlfipV 

- iio^xe v 


m h — 


lou^xe v 






V ev ) 


RT3 ( h _^ hh\ 
UIll ft * Cul 


Jjrv^n ► ^L^ji ) 


Dill ft ^ UC7 1 


VlT)(h — f /-yr 7/™^ 






10 2 


2.74 x 10~ 5 


1.0000 


1.137 x 10~ 5 


0.9999 


<5 = 


3 


10 3 


0.02647 


0.9642 


0.01064 


0.9360 






10 4 


0.7196 


0.02621 


0.1639 


0.01442 






1 


0.59 


0.15 


0.21 


0.10 


<5 = 


2 


10 


0.74 


2.57 x 10~ 3 


0.1662 


1.06 x 10~ 3 






10 2 


0.74 


2.58 x 10~ 5 


0.1663 


1.06 x 10~ 5 



perform the sum over the KK states (c.f. Eq. ([10))) we get for the invisible decay width of 
the Higgs boson, 



> mh 

Tf 105 



m 



u 



(m h RY 



- if 

m,, 



(6 = 2) 



(24) 



(m h R) 3 . (6 = 3) 



The Higgs boson can also decay into the SM decay modes h — > bb and h — > W*W. We 
compute the SM Higgs decay widths using HDECAY Q]. It is well known that for SM 
Higgs boson mass nih < 135 GeV, the decay h — > bb dominates, whereas above 150 GeV, the 
dominant decay mode is h — > WW*. These two decay channels compete with each other in 
the intermediate mass region 135 GeV < < 150 GeV. 

The decay branching ratios of BR(/i — > bb) and BR(/i — > vlVr) as a function of 1/R are 
shown in Fig. Ufa) for = 115 GeV and = 150 GeV for 5 = 2, 3. The exact values of 
the branching ratios are listed in Table [H] Due to the enhancement from the large number 
of KK modes that are accessible at smaller 1/R, the invisible decay mode becomes dominant 
as is the case, for example, for 6 = 3, 1/R ^ 10 3 eV. Fig. Hfb) shows the branching ratios as 
a function of 1/R for 5 = 2,3. 
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FIG. 4: Decay branching ratios as a function of 1/R and rrih for J2i,e \ m u\ 2 = (0-05 eV) 2 . 

Even though the h — > WW* decay mode is comparable to h — > bb for > 150 GeV, 
this mode is harder to be detected, both in the leptonic and hadronic decay modes of the 
WW*, given the suppression in the decay branching ratio in the former and the dominance 
of the background in the latter. Therefore we will not consider the h — > WW* decay mode 
further in this study. 

C. Monte Carlo analysis 

In this section, we perform a Monte Carlo study to show how to detect such a signal 
at the LHC, for Normal, Inverted and Degenerate mass schemes. In order to establish the 
LHC range of sensitivity, we consider two options for the integrated luminosity, 100 fb" 1 and 
1000 fb _1 (SLHC). Since the signal cross section is very small for 5 = 2 (c.f. Fig. we will 
focus on 5 = 3 in the remainder. 

As is evident from Table ITT| the invisible decay mode h vlVr is the dominant one for 



15 



1/R?S 10 3 eV. When the Higgs boson decays into the invisible mode, the produced KK state 
neutrino behaves as a massive, noninteracting, stable particle and thus appears as missing 
energy in the detector. Therefore, we can observe only the leptons plus ]fi T , or hadrons plus 
$ T , for which there isn't enough kinematic information to reconstruct the Higgs boson, and 
we can probe the signal only indirectly through the missing energy distribution. On the 
contrary, in the h — > bb decay channel, we can use the bb pair to reconstruct the Higgs boson 
invariant mass; however, this suffers from the small decay branching ratio into this mode for 
1/R < 10 3 eV. In the following, we describe in succession the h — > bb and h — > vjJ^r decay 
modes. 

We find in agreement with Ref. Q], that for our choice of parameters, if only the SM 
Higgs production process 



is considered, followed by the new physics invisible decay of the Higgs boson (h — > v^v 



the signal cannot be distinguished from the SM backgrounds. However, if the new physics 
z/j^-Higgs production process shown in Fig. is also included, we will show here by 
performing a Monte Carlo analysis that the signal event can be detected in certain regions 
of parameter space. 

1. h — > bb mode 

For the bb decay mode, we have the signal process 



for which the experimental signature is the production of £ + bb (where £ + = e + ,fi + ), 
associated with a large missing transverse momentum carried away by . We do not include 
the r channel due to the additional suppression from the branching ratio of r — » nu, pv. In 
this study we will assume that both the b and b jets in the signal event can be tagged with 
a total efficiency of 50%. 



qq -> W* -> iv)h 




qq' — > W* — > £ + h(^ bb)u 



,(n) 
R 



(25) 
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FIG. 5: Representative Feynman diagrams for the SM background processes 

The intrinsic SM background processes (in addition to the W* h process discussed above) 
for this channel are 

tb : qq' -> W* -> 5t(-> 6W(-> £ + z/)) , (26) 
W66 : qq' -> W*(-> tv)g{-^ 66), (27) 

as shown in Fig. E| In addition to the Wbb and £6 processes, there are reducible SM 
background processes containing misidentified or undetected particles (charged leptons or 
jets) which mimic the signal event. However, they are not expected to pose a serious problem 
after imposing kinematic cuts to veto additional jet (or lepton) activity and requiring the 
event to have a large missing transverse momentum. We will therefore only focus on the 
intrinsic SM backgrounds, Wbb and tb. 

To investigate the potential of the LHC to detect such a signal, we take 1/R = 500 eV 
with 5 = 3 and perform a Monte Carlo simulation for rrih = 115 GeV. To compare the 
relevant background event rates to the signal event rate, we shall assume the integrated 
luminosity of the LHC to be 100fb _1 . Here, we will assume a perfect detector that can 
precisely measure the four-momenta of the final state partons. We require the separation 
in AR = yj (A0) 2 + (At]) 2 between any two observable final state partons (not including 
neutrinos) to be larger than 0.4, where A<f) and Ar] are the separation in azimuthal angle 
and rapidity, respectively. We will require the transverse momentum (px) and the rapidity 
(rj) of b, b and £ + to satisfy the following basic cuts: 

p q T > 15 GeV, \r] q \ < 3.0, 
17 
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FIG. 6: Distribution (Number of events/GeV) of pip, and m T a ^ LHC, in the l + bb $t channel, 
for 5 = 3, rrih = 115 GeV, \/R = 500 eV, after imposing the basic cuts specified in Eq. (|28|). 



p e T > 15 GeV, 



< 2.5, 



$ T > 15 GeV, AR > 0.4. 



(28) 



Fig. shows the signal and background distributions after imposing the basic cuts. We get 
only 8 signal events (in all three neutrino mass schemes), and 2.5 x 10 5 background events. 
Since this mode suffers from such a huge background, the signal events cannot be seen at the 
LHC. Imposing a hard cut on the missing energy or increasing the luminosity (SLHC with 
1000 fir 1 ), would not improve this result much. Therefore, due to the small decay branching 
ratio of h — > bb, we conclude that it would be extremely difficult to use the bb mode to 
directly detect the signal at LHC or SLHC. 
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FIG. 7: SM Background process for invisible decay mode. 
2. — > v lVt& (Invisible mode) 

For the invisible mode, we have the signal process 

q + q'^£ + + u R + H{-^ uuf), 

where I = e, fi, r. For the r channel, the r + further decays either into leptons or hadrons. 
For simplicity we only consider the 7r + z/ decay mode (with BR(r — > nu) = 0.11) of r + to 
detect the signal events. Therefore, the collider signature for the signal process is leptons 
+ $ T (for t = e, /i), and tt + + $ T (for t = r). When I = e, fi we will denote this as the 
P. + Ifi T channel, but when i = t we will denote this as the 7r + ^ T channel. The major SM 
background process is the Drell-Yan charged current process 

qq' — > l + v , 

as shown in Fig. [7[ where £ = e, fi, r. 

We study next the collider signature for the normal, inverted and degenerate mass 
schemes. 

a. Normal Hierarchy: Using the Higgs interaction specified in Eqs. (pH]) and (|15jl . we 
perform a Monte Carlo analysis with 1/R — 500 eV and 1000 eV, to illustrate the nature 
of a possible signal at the LHC (a p — p collider with 14 TeV center-of-mass energy). We 
consider an integrated luminosity of 100 fb -1 , and impose the cuts, called Basic and Second 
cuts, as shown in Table IHIl for lepton+$ T and n + $ T channels. 

We summarize the number of the signal events (S) and background events (B), after 
imposing the Second cuts, in Table IIVI along with the ratio of the number of signal to 
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TABLE III: Kinematic cuts used to select events 









Basic cuts 


Vt 


> 15 GeV 


p T > 15 GeV 




$T 


> 15 GeV 


$ T > 15 GeV 






if 


< 2.5 


\rf\ < 3.0 


Second cuts 


$T 


> 400 GeV 


$ T > 400 GeV 



TABLE IV: Number of signal (S) and background (B) events for the normal hierarchy scheme, 
with 1/R = 500 eV, at the LHC with an integrated luminosity of 100fb _1 . The kinematic cuts 
listed in each row are applied sequentially. 







7T+ +$ T 




t = e 


£ = H 


S 


B 


S/B 


s/Vb 


s 


B 


S/B 


s/Vb 


Basic cuts 


26 


1691 


1717 


1.07 x 10 9 


1.6 x 10~ 6 


0.05 


191 


0.17 x 10 9 


1.12 x 10" 6 


0.015 


Second cuts 


6 


393 


399 


2432 


0.16 


8.1 


58 


27 


2.15 


11.2 



background events (S/B) and the statistical significance of the signal (S/\fB). In Fig. [HI we 
show the missing energy distribution for 1/R = 500 eV and 1000 eV, after imposing both 
the Basic and Second cuts, for £ + ^ T in (a), and n + ]fi T in (b). 

In the £ + ]fi T channel, Fig. IHfa) shows that for 1/R = 500 eV the signal can be 
distinguished from background at large $ T . This is possible because the density of KK 
states (dn in Eq. (fTU|) ) increases with increasing mass of the KK state, due to which the $ T 
distribution tends to peak at higher $ T compared to the background. Even though this 
effect is somewhat reduced due to the rapidly falling PDF's, the signal still tends to peak at 
larger $ T . 

In the tt + $ t channel, as shown in Fig. Efb), the r + momentum in the signal is harder 
compared to the background since it is balanced against the heavy z/^-Higgs system due 
to which the resulting 7r + from the r + decay is also harder. Therefore after imposing the 
hard $ T cut (Second cuts), more background events will be cut away compared to the signal 
events leading to a significance of 11.2 a, as shown in Table HVl 
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FIG. 8: Missing energy distribution (after both Basic and Second cuts) for the normal hierarchy 
scheme, with 1/R = 500 eV and 1000 eV. (a) is for £ + ]£ T channel, and (b) for tt + $J t channel. 

b. Inverted Hierarchy: We follow a similar procedure and perform a Monte Carlo using 
the Higgs interaction specified in Eqs. (fb"f|) and (|16|) . We use the same kinematic cuts given 
in Table UTTl for lepton+^ T and 7r + $ T . The number of signal and background events after 
imposing these cuts are summarized in Table The ratio of the number of signal and 
background events (S/B), as well as the statistical significance of the signal (S/yB), are 
also shown in Table |3 In Fig. we show the missing energy distribution after imposing 
both the Basic and Second cuts, for i + $ T in (a), and w + $ T in (b). As can be seen from 
Table for 1/R = 500 eV, there is a significant deviation from SM background in both the 
l + $ T and n + ]£ T channels. 



We note that for both the normal and inverted mass schemes, the statistical significance 
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TABLE V: Number of signal (S) and background (B) events for the inverted hierarchy scheme, 
with 1/R = 500 eV, at the LHC with an integrated luminosity of 100fb _1 . The kinematic cuts 
listed in each row are applied sequentially. 







7T+ + $ T 




£ = e 


l = [i 


S 


B 


S/B 


s/Vb 


s 


B 


S/B 


s/Vb 


Basic cuts 


3372 


1649 


5021 


1.07 x 10 9 


4.69 x 10~ 6 


0.15 


183 


0.17 x 10 9 


1.08 x 10~ 6 


0.014 


Second cuts 


767 


375 


1142 


2432 


0.47 


23.2 


55 


27 


2.04 


10.58 



TABLE VI: Number of signal (S) and background (B) events for the degenerate mass scheme, with 
1/R = 3000 eV, at the LHC with an integrated luminosity of 100fb _1 . The kinematic cuts listed 
in each row are applied sequentially. 





£ + + $ T 


7T+ + $ T 




£ = e 


l = H 


S 


B 


S/B 


s/Vb 


s 


B 


S/B 


s/Vb 


Basic cuts 


6049 


6065 


12114 


1.07 x 10 9 


1.1 x 10~ 5 


0.37 


672 


0.17 x 10 9 


4.0 x 10~ 6 


0.05 


Second cuts 


1376 


1380 


2756 


2432 


1.13 


55.9 


202 


27 


7.48 


38.87 



for distinguishing the signal from background becomes marginal for 1/R ■> 1000 eV, due to an 
overall suppression of the signal rate in the entire $ T range, as suggested by Figs. (jH)) and (jHJ). 

c. Degenerate: As already pointed out, we expect the constraints on 1/R from neutrino 
oscillation and unitarity to be stronger in the degenerate mass scheme, and we therefore 
choose to present results for 1/R = 3000 eV, as this yields the same order of magnitude 
total cross section as in the normal and inverted mass schemes. We perform a Monte Carlo 
study with the Higgs interaction given by Eqs. (fT^j) and (|17|). The numerical results are 
summarized in Table. IVII and Fig. and in the latter we also show the distributions for 
1/R = 6000 eV. In the degenerate mass scheme, there are about the same number of signal 
events in the e + , /i + and r + channels. 
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FIG. 9: Missing energy distribution (after both Basic and Second cuts) for the inverted hierarchy 
scheme, with 1/R = 500 eV and 1000 eV. (a) is for £ + ]£ T channel, and (b) for ir + ]ft T channel. 

D. Discovery potential 



In order to study the potential of the LHC to distinguish the signal from background 
events, we compute the Significance 

S O signal^ 



s 



B 



v^B J a bkgd C 



(29) 



where £ is the integrated luminosity. To unambiguously establish the new physics scenario 
we are considering, it is desirable to have Sb > 5. For each mass scheme, we plot in Fig. ^2 
Sb as a function of 1/R, at the LHC for £ = 100 and 1000 fb _1 . In this, we have imposed 
both Basic and Second cuts listed in Table. ITTT1 for the l + $ T and n + $ T channels. 
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FIG. 10: Missing energy distribution (after both Basic and Second cuts) for the degenerate mass 
scheme, with 1/R = 3000 eV and 6000 eV. (a) is for £ + $J T channel, and (b) for 7r + $y channel. 

We infer from Fig. ^2 that the C. + $ T and ir + $ T channels lead to a similar significance, 
with l + $ T being somewhat better in the case of the inverted mass scheme. If 1 /i? is about 
900 eV (for normal or inverted) or about 8 KeV (for degenerate) we may expect Sb in the 
2 — 5 a range at the LHC Here, it is worth pointing out that if we strictly impose the unitarity 
and oscillation constraints on 1/R that we derived in Ref. we would expect a poor Sb 
at the LHC with C = 100 fb^ 1 . However, we had pointed out that owing to the dependence 
on the cutoff for 5 = 3, those constraints are uncertain to some extent. Therefore, in the 
event that the constraint on 1/R is relaxed somewhat, we may expect to see a signal at the 
LHC. 



The signal and background event rate estimates discussed thus far are subject to PDF 
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FIG. 11: Significance (after Basic and Second cuts) as a function of 1/R (5 = 3) for the normal, 
inverted and degenerate mass schemes at LHC with an integrated luminosity C = 100 and 1000 fb -1 . 



uncertainties 
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21] . Since both our signal and background processes are mediated through 



a virtual W* production, we can estimate this uncertainty by scanning over the 41 sets 



of CTEQ6.1 parton distribution functions 



2l| to calculate the next-to-leading order W* 



production cross section. We find that the uncertainty in W* (Drell-Yan) production, over 
the Q 2 region of interest, is about 2 % at the LHC energy. Since both our signal and 
background processes go through a W*, this uncertainty largely cancels and yields about 
the same S/B value. Furthermore, when data becomes available at the LHC, the PDF 
uncertainties will be further reduced from global analysis. Another way is to use the side- 
band method in order to estimate the SM background rate in the signal region, enabling us 
to minimize the PDF uncertainties. For example, in our case we could use the background 
dominated $ T < 400 GeV region as the side-band to estimate the SM background rate. 
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E. Distinguishing mass schemes 



We pointed out in Sec. Ill Cl that the neutrino oscillation data presently leave undetermined 
the sign of Am 2 tm , in addition to the absolute scale of neutrino masses. Measuring the sign 
of Am 2 tm determines directly whether the normal or the inverted hierarchy is realized in 
nature. While presently no experiment is capable of achieving this, the proposed long baseline 
neutrino experiments could potentially determine this sign by making use of matter effects 



in the earth |22j. Here, it is interesting to ask if we can instead use the collider observables 
that we have been considering in order to distinguish between the normal, inverted and 
the degenerate mass schemes, given an excess above SM background at the LHC that is 
consistent with bulk right handed neutrinos as is being discussed here. We define M{fi + $ T ) 
and A/"(e + $ T ) to be the number of (/i + $ T ) and (e + $ T ) signal events, respectively, after 
the Second cut. As a first step, we assume that N(n + $ T ) and M(e + $ T ) can be determined 
by subtracting out the estimated number of background (p + lfcp) and (e + $ T ) events, from 
the actual number of events seen in the collider. We find a suitable discriminant that could 
potentially distinguish between the neutrino mass schemes to be 

_ N^ + $ T )-N{e + $ T ) 



It can be shown from Eqs. (|T|). ()14j) and (jl9j) that 

{ifR) 
£ 

(IfR) 

where the upper (lower) sign is for the normal (inverted) hierarchy.^ Therefore, 



Af(e + $ T ) oc [ml + (0.531) 2 AmL ar 

A% + $ T ) oc 7T ^ [ml + (0.599) 2 AmL ar ± (0.707) 2 Am 2 J , (31) 



_ ±0.5Am 2 atm 

and for normal (inverted) hierarchy A^ e > (A^e < 0). 

In Fig. ^1 we plot A^ e as a function of m\. We note here that for the inverted mass 
scheme, the smallest value that m x can take is 0.05 eV, as can be seen from Fig. [TJ We 



■ For the Second cuts with i5 = 3, we find, say from Table IVT1 the constant of proportionality in Eq. to 
be 3.726 x 10 11 eV fb. 
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FIG. 12: Aue as a function of mi. w4^ e > (A^ e < 0) implies normal (inverted) hierarchy and 
•A- ne ~ implies degenerate mass scheme. The error-bars are for the LHC with C = 100 fb _1 and 
l/R = 500 eV. 

see from Fig. El that it is indeed possible to distinguish between the normal and inverted 
hierarchies. The vertical error-bar shows <L4 Me , the statistical accuracy with which A^e can 
be measured at the LHC, for a luminosity C = 100 fb" 1 . We estimate SA^e by calculating 
the number of e and \i events using Eq. (|3ip. For m,\ small, we see that there is excellent 
discriminating power of many standard deviations to determine whether the normal or the 
inverted mass hierarchy is realized. However, as mi increases the number of e, fi and r 
events become approximately equal and the discriminating power diminishes. 

The absolute scale of neutrino masses (say mi) can also be determined by measuring 
A^e, but only if mf is not too large compared to Am 2 . This is because, a given 
(statistical) uncertainty SA^ e corresponds to an uncertainty on the mass given by 8m\ = 
(0.25Am2 tm )(&4 /je /.Aj e ). Thus, for m\ > Am 2 the masses become approximately 
degenerate, M{e + $ T ) ps A/"(/x + $ T ) ~ M{r + # T ), due to which A^ e w 0, and the 
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uncertainty 8m\ becomes worse. From Fig.^Jwe see that we can determine mi by measuring 
A^e, and the horizontal error-bars shows the uncertainty in the inferred value of mi due to the 
(statistical) uncertainty SA^e on the measured A^ e . This might provide a means to measure 
the absolute scale of neutrino mass (mi) which cannot be done through neutrino oscillation 
experiments. For example, in the normal mass scheme if mi = 0.05 eV, at the LHC we 
would measure an asymmetry A^e = 0.2 ±0.03, where the error is the statistical uncertainty 
for £ = 100 fb _1 . From this measurement, we infer from Eq. (|3*2*|) that mi = 0.05lo!oo5 e V, 
and this error is shown in Fig.^Jas the horizontal error-bar. For mi > 0.2 eV the horizontal 
error-bars get progressively larger and we do not show them in the figure to reduce clutter. 
If mi is indeed large in nature, it follows from Eq. (J22)) that from the measured bulk neutrino 
production cross section, we can at best only determine the ratio m 2 /(1/R) 5 using Eq. (|19|) . 

It is interesting to note that the neutrino oscillation probability from an active species to 
the (sterile) heavier KK modes is also proportional to m 2 /(l/R) s [ljj. Therefore, it appears 
that in the degenerate mass case, it would not be possible to disentangle m and 1/R even 
if we observed the collider signature that we have been considering; this would still be the 
case if we also observed a finite oscillation probability into the heavier KK states. 

In Eq. (jl2j) . for simplicity, we assumed Z e3 to be zero. If we had not done so, one might 
wonder if A^e has any sensitivity to the small angle associated with / e3 . We do not expect 
this to be measurable since the change in the numbers of e and \i events would be much 
smaller than the statistical uncertainties involved. It is also worth mentioning that the r-e 
(or T-fi) asymmetry does not add any information in probing the neutrino masses. 



IV. CONCLUSIONS 



We consider a theory with right-handed neutrinos propagating in 5 large extra dimensions, 
with radius R, and coupled to the left-handed neutrino by the Higgs boson Yukawa coupling. 
Such a theory naturally explains the smallness of the neutrino mass scale, in addition 
to addressing the gauge hierarchy problem in the standard model (SM). Presently the 
combination of neutrino oscillation data leave undetermined the absolute neutrino mass 
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scale m v which lead us to present our results for the normal, inverted and degenerate mass 
schemes. In addition, the observables that we analyze depend on 8 and 1/R of the extra- 
dimensional theory. We take the fundamental scale of gravity M* to be the electroweak scale 
M EW ~ 10 3 GeV. 

We showed in our previous work that neutrino oscillation and unitarity in Higgs-Higgs 
scattering places a lower bound on 1/R, particularly strong for 5 > 3. We thus restricted 
our consideration to theories with 5 < 3, and we pointed out that the bounds we derived 
were somewhat uncertain for 5 = 2 and 3, since the observables depended on the cutoff (M*) 
of the theory. Here, we present our results choosing 1/R to be around this lower bound, 
although not imposing the bound strictly due to the uncertainty just mentioned. 

We show that the Higgs boson production at a collider and its decay can be enhanced 
significantly owing to the large number of Kaluza-Klein (KK) neutrino states that can be 
produced in the final state. To probe this theory we consider the signal process 

qq' ^W* ^£+h4 in) 

at hadron colliders, such as the Tevatron and the LHC The production cross section is 
shown in Fig. (jHJ) for 5 = 2, 3. We find for instance, for S = 3, 1/R ~ 1 KeV, we could detect 
a signal at the LHC with a luminosity C = 100 fb~ l . 

In addition to the SM decay mode of the Higgs boson h — > bb, we also consider the new 
invisible decay mode (h — > ^l^r)- As shown in Fig. (jl]), the invisible decay mode dominates 
for smaller 1/R. 

We perform a Monte Carlo analysis of the signal process along with the SM backgrounds, 
and find that the Higgs boson decaying invisibly, leading to a £ + $ T (£ = e, \i, r) signature, is 
more promising compared to the h —>■ bb decay mode. We apply the cuts shown in Table HTD 
to enhance the signal relative to the background, and show the resulting significance in 
Fig. We find a significance in the 2 — 5 a range for 5 = 3 and 1/R ~ 900 eV for 

normal and inverted mass hierarchies, and 1/R ~ 8000 eV for the degenerate mass scheme. 
We only consider the positive charged lepton and the 7r + v decay mode of the r + . Including 
the negative charged lepton and the pv decay mode of the r will improve the situation 
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somewhat. 

Finally, we point out that if a positive signal is found that is compatible with the extra- 
dimensional hypothesis considered here, the asymmetry in the number of \i versus e events 
can be used to distinguish between the neutrino mass schemes and to determine the absolute 
neutrino mass scale, and the accuracies with which these can be determined are shown in 
Fig. (|12|). The collider might be a unique place to determine this information as it is not 
available from neutrino oscillation experiments. 
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